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We discuss the dynamical behavior of the entanglement between the internal and the external de- 
grees of freedom of a trapped atom in electromagnetically-induced transparency (EIT) laser cooling. 
It is shown that essential features of the intricate entanglement dynamics observed in full numerical 
simulations of the underlying quantum master equation can be understood in terms of a two-state 
model on the basis of Landau-Zener splittings in the atom-laser field Hamiltonian. An extension 
of this model to an effective non-Hermitian Hamiltonian is constructed which describes the decay 
of entanglement by spontaneous emission processes. We also discuss schemes for the control of 
entanglement and demonstrate that a permanent entanglement can be imprinted on trapped atoms 
through a rapid switch off of the driving fields. Finally, we point out fundamental distinctions be- 
tween the entanglement created in EIT cooling and in the cooling scheme based on velocity-selective 
coherent population trapping. 
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I. INTRODUCTION 

Entanglement is a characteristic property of the states 
of composite quantum systems and represents a key fea- 
ture for many proposed applications of quantum mechan- 
ics. Since quantum systems are often strongly affected 
by decoherence through the interaction with an environ- 
ment pQ, the dynamics of entanglement in the presence 
of dissipative processes plays an important role in many 
contexts [2]. One of the challenges is, for example, to find 
ways to protect the coherence of quantum states against 
external noise and to maintain entanglement for an ex- 
tended period of time. 

Building on the ideas presented in [3] we investigate 
here in detail the dynamics of the entanglement that 
arises between the internal electronic and the external 
translational degrees of freedom of trapped atoms in laser 
cooling. The cooling scheme investigated employs the 
cooling of the vibrational motion of the atom in an ex- 
ternal potential by means of electromagnetically-induced 
transparency (EIT) using two counter-propagating and 
suitable tailored laser fields. In EIT cooling entangle- 
ment is a fundamental requisite for cooling, but fades 
once cooling into the vibrational ground state is achieved. 
The importance of entanglement for the cooling process 
permits to imprint permanent entanglement between the 
atom's internal and external degrees of freedom. 

The situation studied here is reminiscent of the cool- 
ing mechanism based on velocity-selective coherent pop- 
ulation trapping (VSCPT) which has been considered by 
Arimondo fj] who studied the entanglement between mo- 
mentum eigenstates of the free atom and its internal elec- 
tronic degrees of freedom. By contrast, here we study 



the entanglement between vibrational eigenstates in the 
external potential and the atom's internal electronic de- 
grees of freedom. In the ultimate cooling limit, the de- 
gree of entanglement reached is very different from that 
in VSCPT. In the latter case a maximally entangled state 
appears when the translational motion is fully cooled [3] . 

The degree of entanglement in the atomic density ma- 
trix p(t) will be quantified here by means of the negativity 
[5] which is defined as 



AT(i) = AT(p(£)) = |(||p T (t)|| - 1 



(1) 



'Electronic address: maryam.roghani@physik.uni-hamburg.de 
^Electronic address: breuerSSphysik.uni-freiburg.de 
JElectronic address: helm@uni-freiburg.de 



Here, p T (t) represents the partial transpose of p(t) with 
respect to either the internal or the external degrees of 
freedom, and ||^4|| = Tr|A| denotes the trace norm of 
an operator A. The negativity N{p) is a nonnegative, 
convex function and represents a useful measure for the 
amount of entanglement of quantum states which is easy 
to compute. It represents an entanglement monotone, 
vanishes for separable states and takes on a maximum 
value if p is a pure, maximally entangled state. Per- 
forming extensive numerical simulations of the underly- 
ing quantum master equation that describes the dynam- 
ics of the atomic density matrix, we are able to study the 
complex time evolution of the negativity in EIT cooling. 

The paper is organized as follows. After a brief intro- 
duction of the EIT-cooling scheme and of the quantum 
master equation used to model this scheme, we present 
in Sec.|n]the intricate dynamical behavior of the negativ- 
ity which is observed during the vibrational cooling pro- 
cess [2) and show that the basic features of the entangle- 
ment dynamics can be traced to Landau-Zener splittings 
(avoided crossings) [SHE] in the spectrum of the atom- 
laser field Hamiltonian. We further construct an effec- 
tive Hamiltonian that describes the behavior of entangle- 
ment under dissipation induced by spontaneous emission 
processes. In Sec. Ill we present results of numerical 



simulations of the full master equation of EIT cooling, 
and compare these with the predictions derived from the 



2 



Landau-Zener model. Furthermore, we discuss here the 
control of entanglement through sudden and adiabatic 
switching off the driving laser fields. Finally, we draw 
our conclusions in Sec. HVI and discuss in some detail the 
differences in the entanglement dynamics of the cooling 
schemes of EIT and VSCPT. 



II. THEORY 




A. Theoretical model and master equation 



FIG. 1: Level scheme of the trapped atom. 



In a trapped atom system the Lamb-Dicke parameter 
i] = kj \J1mui ih.=Y) controls the probability of changing 
the vibrational quantum number n in electronic transi- 
tions [9]. For small r/ the photon-recoil energy k 2 /2m 
is an inconsiderable fraction of the trap vibrational fre- 
quency to. In the Lamb-Dicke regime, r\\/n + 1 <C I, 
electronic transitions with a change of the vibrational 
level are active, primarily in the first sideband transi- 
tions, An = ±1, their probability being proportional to 
(?i + f)?7 2 , where n is the vibrational quantum number. 
A rapid removal of vibrational energy of a trapped atom 
occurs for suitable experimental parameters [TOUT!?] , i.e., 
a difference of ac-Stark shifts of the two ground states 
of the order of uj. Under these conditions, with blue de- 
tuning of two counter-propagating lasers, a vibrationally 
excited atom can be driven into states with mean values 
(n) near zero, as has been demonstrated in experiments 
with atomic ions |13) . 

We have explored this concept in a model based on the 
full quantum master equation 



df 



p{t) = -\[H\p{t)]+Cp{t) 



(2) 



for the interaction picture density matrix p(t) represent- 
ing the combined state of vibrational and electronic de- 
grees of freedom of the trapped atom. The total Hamil- 
tonian 



H* 1 = H r 



(3) 



of the model consists of three parts: H cm = uia^a de- 
scribes the vibrational degree of freedom of a harmoni- 
cally trapped atom with raising and lowering operators 
<v and a, and H c \ = A(|l)(l| + |2)(2|) represents the elec- 
tronic degree of freedom with detuning A. The electronic 
states are denoted by \i), i — 1,2,3, and form a A-type 
level structure sketched in Fig. [T] The Hamiltonian 



^t = |e ifex |3)(l| + |e- 



|3)(2|+h.c. (4) 



models the interaction with two counter-propagating 
laser beams, the most efficient configuration for cooling, 
where x denotes the atomic position and k the wavenum- 
ber (hi ~ k 2 — k). Note that H^ nt depends on the Lamb- 
Dicke parameter r\ through kx = r\(a) + a). 



For EIT cooling the detunings are chosen positive 



Ai = U3 L i - U3 31 > 0, 



A, 



U L 2 ~ ^32 > 0, (5) 



and equal, Ai = A2 = A. The latter condition implies 
that in the dressed state frame the non-interacting states 
|1) and 1 2) are at equal total energy and separated from 
the excited state |3) by the detuning A. 

Finally, dissipation caused by spontaneous emissions 
is modeled by the Lindblad superoperator (see, e.g., 
Ref. P) 



3=1,2 q=± 



a-jqpa 



(6) 



where T denotes the spontaneous emission rate and 
0j+ = b')(3|e lfcr , o-j- = \j}(3\e~ lkx . The Lindblad oper- 
ators Uj + describe the transitions from the excited elec- 
tronic state 1 3) to the ground states \j = 1, 2) with a pho- 
ton emitted parallel to the harmonic oscillator axis, while 
the <jj- describe the transitions from |3) to \j = 1,2) 
with a photon emitted counter-propagating to the har- 
monic oscillator axis. This is a good compromise short 
of doing a full 3D simulation as the recoil of the transla- 
tional motion is zero for emissions perpendicular to the 
trap axis, while it leads to Doppler heating or cooling for 
emissions along the trap axis. 

Typical results obtained by numerically solving Eq. ^ 
[31 Q2] are given in Fig. [2^, which gives the dependence 
of the mean vibrational quantum number as a function 
of time for 77 = 0.1, for three values of the excited state 
width r. The analysis of the radiation field which is 
scattered from the atom during the process of cooling 
[T2"] showed that it is the interplay of stimulated Ra- 
man transitions between vibronic states and spontaneous 
emissions which controls the speed and degree of cooling. 

As the Raman transitions are active in a space of in- 
separable wavefunctions it appeared useful to us to also 
explore the temporal development of entanglement dur- 
ing the cooling process. As apparent from Fig. [2]d an in- 
tricate behavior of the negativity is observed during such 
a vibrational cooling process. We see a rapid build-up of 
entanglement, a subsequent exponentially damped oscil- 
lation, followed by a precipitous drop and a power law 
behavior and, finally, a rebirth of entanglement towards 
the non-equilibrium stationary state. 
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FIG. 2: a) Dynamics of vibrational cooling for the initial vi- 
brational level rii — 10 for rj — 0.1 and for three values of the 
excited state width V. b) Dynamics of the negativity during 
the cooling process. The top scale gives the mean vibrational 
quantum number for the 6 MHz result. The dotted line rep- 
resents a power law decay, M oc t~ 3 / 2 . The parameters are 
gi = 1.34, pi/32 = 10, oj = 0.03 and A = 15 in units of 2tt MHz. 



We show below that this complex behavior is related to 
Landau-Zener splittings or avoided crossings in the spec- 
trum of the atom-laser field Hamiltonian. The dominant 
source of entanglement is a near degeneracy of vibronic 
states, which appears as a consequence of the ac-Stark 
shift, suitably chosen for EIT-cooling. Two-photon tran- 
sitions associated with a change in vibrational quantum 
number occur and lead to a periodic oscillation of entan- 
glement, damped by spontaneous emission. 



B. Landau-Zener model for entanglement dynamics 

1. Near degeneracies of vibronic states 

To introduce the Landau-Zener model we first neglect 
spontaneous emission and set the Lamb-Dicke parameter 
to zero. Then the interaction Hamiltonian is 



HL = Y |3)<l| + f 13} <2 



.92 



(7) 



EIT-cooling is most efficient with regard to the cooling 
limit if one of the lasers is very much stronger than the 



other one, such that the difference of ac-Stark shifts expe- 
rienced by the ground states is approximately equal to to. 
This situation is illustrated in Fig. [3j which gives the 
dressed states of the harmonic oscillator, spatially sepa- 
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FIG. 3: The dressed atomic states in the absence of inter- 
action for gi 3> t/2 and an ac-Stark shift of magnitude uj, dis- 
placed along the horizontal axis for clarity. 

rated for clarity. The assumption is made that laser one 
is very much stronger than laser two, thus the ac-Stark 
shift between |1) and |3) repels these two states while the 
shift of state |2) is negligible. For trapped neutral atoms 
we typically have A ^> w, the separation between the 
dressed ground states and the excited state being much 
larger than it appears from Fig. [3] 

In the dressed state picture, the eigenstates of the 
Hamiltonian H c \ 



H? lt which are given by 



(H el + H? nt )\<f> i )=s i \<i >i ) 



have the energy eigenvalues 

1 

e 2 = +A 



n+A 



(8) 

(9) 
(10) 

(11) 



where we have abbreviated = \J g\ + g\ + A 2 . The 
corresponding eigenvectors \<pi) are given by 



1 

yj2VL (Q - A) 

1 (- 



.91 



y/2Q. (VI + A) 



( ffl |l)+ 52 |2) + (0-A)|3)X,L2) 

.92|l>+.9i|2)), (13) 
( 5l |l)+ 92 |2)-(0 + A)|3)X.14) 



With the choice gi^>§2 we have |</>2) ~ |2), while the state 
carries a small contribution of the excited electronic 
state 1 3). 

For vanishing Lamb-Dicke parameter, no coupling be- 
tween electronic states exists via sideband transitions. In 
this case the three electronic states are always tied to a 
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FIG. 4: Energy eigenvalues as a function of the parameter A. 
The lowest Landau-Zener crossing appears for r\ = at A = 1. 
Energies are given as Ei : „ — A. Left: n = 0, right: interacting 
states for a Lamb-Dicke parameter n = 0.1 . 



specific vibrational level. The eigenstates of the total 
system described by the Hamiltonian 



(15) 



are given by tensor products of the dressed electronic 
states \<f>i) and the vibrational eigenstates |n), 



\4>i,n) = \4>i) <g> \n) 



(16) 



with the corresponding eigenvalues £i.„ = E{ + nuj. When 
the ac-Stark shift is adjusted such that the two-photon 
detuning between Stark-shifted states equals the trap fre- 
quency, 



(17) 



then the states \<f>i,n— 1) with energy £i !Tl _i and \(f>2,n} 
with energy e 2 n become degenerate, as is illustrated in 

Fig.m 

For a more general description we consider Rabi fre- 
quencies gi(X) which vary linearly with the parameter A, 
51(A) = Xgi and 32(A) = A 32- Under the resonance con- 
dition ( 17) the energies of the pairs of states \4>i, n— 1) and 
\4>2, n) then cross at A = 1, as is shown in the energy di- 
agram of Fig. [4j It is known from the Landau-Zener sce- 
nario that this exact degeneracy is lifted under a generic 
perturbation. In fact, for a nonvanishing Lamb-Dicke 
parameter, rj > 0, the levels exhibit an avoided crossing 
(Landau-Zener splitting) with a certain energy gap AE, 
as is shown in Fig. [4]b. 

To estimate the energy gap AE we consider the first 
order approximation of the interaction Hamiltonian H? nt 
which is given by 



the unperturbed eigenstates \(f>i,n — l) and \4>2,n) one 
obtains for the energy gap 



where 



AE = E+ -E_ = 2\A\, 



A={4> u n-l\H^\<t>i,n). 



(19) 



(20) 



Under the condition A 3> g\ , g 2 this matrix element may 
be approximated by 



A « —vt]\fn 



9i92 
2A ' 



(21) 



With an appropriate choice for the phase of the unper- 
turbed states the eigenstates |x±) can be written as 

\x±) = ^{\4>un-l)±\^,n)). (22) 



2. Build-up of entanglement 

The key features of the dynamical generation of entan- 
glement between the electronic and vibrational degrees 
of freedom can be understood with the help of a simple 
approximation based on the Landau-Zener picture intro- 
duced above. The general solution of the Schrodinger 
equation 



-mom)) 



(23) 



in the space spanned by the Landau-Zener pair of states 
\X±) is given by 



|*(t)) = a+ e- 1 ^| X+ >+a_e-^- t | X _). 



(24) 



For the initial condition \^(t = 0)) = \(f>2, n), for example, 
we have a± — ±l/y/2 which leads to 



|tf(t)) = -[e- [E+t -e- iE -* 



E+t 



\(f>i,n-l) 
B - t ]|^ 2 ,n). (25) 



The right-hand side represents the Schmidt decomposi- 
tion of the wave function with Schmidt coefficients 



Oil 
«2 
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\E + t 
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\E + t 
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iB_ 1 1 



(26) 



"7, 



ff i|3)(l|- 5 i|l)<3|- 52 |3)(2|- 



ff 2 |2)(3| 
(18) 

The perturbed eigenstates right at the avoided crossing 
(A = 1) are denoted by |x±) with corresponding eigenen- 
ergies E± . Employing degenerate perturbation theory to 



The negativity of the state ( 25 ) is given by the product 



of the Schmidt coefficients |14) and, hence, we obtain for 
the negativity of an isolated Landau-Zener state pair 

Afuz(t) = aia 2 =-\wn.(AEt)\. (27) 



5 



Of course, the same result is obtained for the initial state 
\^(t = 0)) = |0i,n—l). Thus, we see that within the two- 
state Landau-Zener approximation the entanglement dy- 
namics at the avoided crossing is determined completely 
by the spectral gap AE, leading to a periodic oscilla- 
tion of entanglement between zero and \ with the period 
?cnt = tt/AE. It should be noted that considering only 
a single isolated state pair is a valid approximation only 
as long as the splitting AE is much smaller than the 
separation lj between different Landau-Zener state pairs. 
According to Eqs. (191 and ( 21 1 this condition sets an up- 



per bound for the Lamb-Dicke parameter 77 and for the 
vibrational quantum number n when considering multi- 
ple Landau-Zener state pairs. 



3. Decay of entanglement 

Up to now we have discussed the case T = 0. There is 
a small contribution by the excited electronic state |3) in 
the dressed state 1 0i) which according to Eq. (14 1 is of 
magnitude yJojJVL 



y/uj/A and opens the decay channel 
due to spontaneous emissions. In order to account for the 
effect of spontaneous emissions within the Landau-Zener 
model we consider the Hamiltonian 



H° eS = H°-i-\3){3\. 



(28) 



This is the non-Hermitian effective Hamiltonian corre- 
sponding to the master equation ([2| with the Lindblad 
dissipator |6]) which describes the coherent dynamics of 
the Hamiltonian H° including the finite decay widths 
of the dressed states caused by spontaneous emission 
processes. By use of this effective Hamiltonian we ne- 
glect real emission processes which lead to transitions 
into other Landau-Zener state pairs. Thus, we neglect 
here the gain of populations of other state pairs which 
are assumed to yield only negligible contributions to the 
negativity. We will see later in our numerical simulations 
of the full master equation that this assumption leads to 
a good approximation for the entanglement dynamics in 
the early phase of the cooling process. 

The complex eigenvalues of H® s are given by Ei >n — 
£i + noj, where 

ei = ir/2) + y(A + ir/2) 2 + g\ + 9I 

£2 = A, 

S3 = ^(A-iT/2)-^/(A+ir/2j^ 



9\ 



9l 



The corresponding eigenstates are of the tensor product 
form \<t>i,n) = \<f>i) <g> \n), where 



= JVi( fl i|l)+ fla |2)+2(ei-A)|3)) 

= AT 2 (-<7 2 |l) + <?i|2)), 

= AT 3 (<?i|l)+<72|2)+2( £3 -A)|3) 



with the normalization factors 



Ni 



92 



■4k,,- Al 2 



(29) 



The imaginary parts of the complex eigenvalues £j de- 
scribe the widths of the respective eigenstates. No width 
appears for the state \<j>2)- Under the condition A » 
9i , 92 one finds 



A 



91 



9l 



.71 



4 A 2 + (r/2) 2 



with the width of the state |0i) given by 



7i 



9i 



9l 



4 A 2 



(r/2) 2 



(30) 



(31) 



We note that the width 71 is much smaller than the width 
of the excited state \<$>z) which is approximately equal to 
L. The small decay path which is active here is marked 
by the red wiggly lines in Fig. [5] 
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02,0) 




first-order coupling 

FIG. 5: Elements which connect to |02, 0) in first order of n. 
The dressed states are shown for A = 1. Two-photon tran- 
sitions in first sidebands are marked as blue arrows, their 
strength given in units of ~ir/gig2/(2A), see Eq. ( |2 1 \ . Vi- 
brationally diagonal transitions for spontaneous decay (red 
wiggly lines) occur at the rate 71 given in Eq. (31 1. 



The dominant width is carried by the eigenstate \4>3). 
This state is well separated in energy from the dressed 
states |0i ) and |0a)- Thus, we can derive again an 
approximation for the negativity by use of degenerate 
perturbation theory involving the states \<f>i,n— 1) and 
|02j ft) • A straightforward calculation analogous to the 



one carried out in Sec. II B 2 then yields the following ex- 
pression for the negativity in the presence of spontaneous 
emission, 



M.z(t) = -g-e" 



7l*/2 



sin vt 



_7i 
2v 



(1 — cosvt) 



(32) 



6 



where v = yj AE 2 — 7^/4. According to this formula the 
entanglement dynamics at an avoided crossing is deter- 
mined by the spectral gap AE and by the width of 
the dressed state \<pi,n). Note that in the limit T — > 
which implies 71 — > the expression ( 32 ) reduces to the 



AMt) 



expression of Eq. ( 27 I 



III. NUMERICAL SIMULATION RESULTS 

In this section we present results of numerical solutions 
of the full master equation ([2| and discuss and interpret 
these solutions on the basis of the Landau-Zener model 
developed in Sec. |II B| We shall demonstrate that the 
Landau-Zener picture is in fact very useful in explaining 
and understanding the essential features of the complex 
entanglement dynamics observed in Fig. [2]d. In this fig- 
ure we see a rapid build-up of entanglement, a subsequent 
exponentially damped oscillation, followed by a precipi- 
tous drop and a power law behavior and, finally, a rebirth 
of entanglement towards the non-equilibrium stationary 
state. In the following we discuss these features and es- 
tablish a connection to the Landau-Zener model. 



A. Build-up of entanglement 

When EIT cooling begins from a state \<f>2,n- m ) with 
n- ln the initial vibrational quantum number, the two-state 
model predicts an initial growth of negativity from zero 
and subsequent oscillations according to Eq. ( p7| and 
( 32 1 . This is borne out by the full simulation and par- 
ticularly apparent for low values of F, see Fig. [2j The 
full simulation also predicts this oscillating negativity 
when r = and in Fig. [6^, and [6]d we compare pre- 
dictions of Eq. (27) with the results of such full simu- 
lations. 



The two-photon Raman transitions which con- 
nect the horizontal Landau-Zener pairs in Fig. [5] are the 
origin for this oscillation. To a much weaker extent 
the two-photon Raman transitions also connect to the 
next-neighbor Landau-Zener pairs, above and below, see 
Fig.[5j This connection is not part of the two-state model, 
thus the small discrepancies between the full simulation 
and the two-state model in Fig. [6] are not surprising. 

Figure [6^ and b give results for two different initial 
vibrational levels, the different periodicity reflecting the 
71-dependence in Eq. (21). The negativity oscillation pe- 
riod is T ont = n/AE. Since AE cx y/n^, the choice of 
different n in values leads to different oscillation periods. 
Note that the full simulation predicts negativities which 
slightly exceed the value of 1/2, the maximal value pos- 
sible for a tensor product space in which one factor is 
two-dimensional. The reason is that the space describ- 
ing the internal degree of freedom of the atom is three- 
dimensional here, the maximal value for the negativity 
being thus equal to 1. However under EIT cooling condi- 
tions the excited electronic state is extremely rarely oc- 
cupied and the actual negativity therefore stays far away 
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FIG. 6: The red lines represent the negativity based on the 
full numerical simulation, 77 = 0.1. Gray lines are the predic- 
tion ( 27 1 based on the Landau-Zener model, using AS = 0.012 



and 0.073 respectively. Parameters are in units of 2ir MHz: 
w = 0.03, gi = 10ff 2 = 1.34, A=15, T = 0. 



from this maximal value. 



B. Decay of entanglement 

When taking into account the small contamination of 
\4>i) by excited state character and its decay width 71, 
we expect decoherence to occur, its magnitude being 
weighted by the amplitude of \4>i) in | In the two- 
state model the trapped atom is initially in |66 2 ) occupy- 
ing the vibrational quantum number |n; n ) with M = 0. 
When the lasers are turned on with A = 1, the atom os- 
cillates between states |</>2, rii n ) and \<f>\, n ln — 1), the nega- 
tivity climbing periodically to a value 1/2. When during 
this cycle the atom finds itself in the state |(/>i, n- m — 1), 
it may suffer relaxation by spontaneous emission to each 
state of the lower Landau-Zener pair, \<f>2 1 n; n — 1) and 
l^ij^in - 2). Both contain a fraction of the electronic 
ground states |1) and |2), but in the Lamb-Dicke limit 
the spontaneous transition to l^^in - 1); indicated by 
the red wiggly arrows in Fig. [5] is stronger by a factor 
oc rf~ 2 . This loss is accounted for within the Landau- 



Zener model by the calculations presented in Sec. II B 3 
In Fig. [7] we compare these predictions of the two state 
model with results from a full simulation at various val- 
ues of r. In this comparison we have used the quantities 
AE and 71 as fit parameters. We see that the behavior 

TABLE I: Fitted parameters used in Fig. [7] in comparison to 



the predictions from Eq. (19 1 and (32) 



r 


AS (fit) 


AE 


7i (fit) 


7i 


100 kHz 


0.0126 


0.015 


0.0016 


0.00063 


200 kHz 


0.0126 


0.015 


0.0016 


0.00126 


600 kHz 


0.0146 


0.015 


0.0062 


0.00377 


6 MHz 


0.0316 


0.015 


0.1141 


0.036 



observed for the negativity mimics that of a damped os- 
cillator, showing a transition from underdamped to over- 
damped motion. As is illustrated in Fig. [7] the simple 
formula ( 32 ) provides an excellent approximation of the 
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FIG. 7: Transition from underdamped to overdamped oscilla- 
tory behavior of the negativity for various values of F, 77 = 0.1, 
ri]„ = 3. Fits of Eq. (32 1 are shown by the gray lines (Table|T|. 



Red lines show the negativity from the full simulation. 



I .0- 



I Pinjm I 



10 12 14 16 18 




n 



entanglement dynamics for 71/2 < AE n , and even for 
intermediate values of the spontaneous emission rate the 
qualitative behavior of the negativity is reproduced. It 
is remarkable that Eq. ( 32 ) correctly describes the tran- 



sition from the underdamped to the overdamped motion 
of the negativity at 7i/2 = AE, compare Figs. [7b and 
[7|i. Using the approximations AE « ^s/ngig^jA and 
71 w rgf/(4A 2 ) this transition is predicted to occur at 
the value T — 8r/y/nAg2/gi, which nicely fits to the re- 
sults of our numerics. 

The interpretation of the damping of negativity is that 
after a spontaneous event the system has forgotten the 
history of coherent oscillation in the former Landau- 
Zener pair and starts to rebuild negativity in the new 
Landau-Zener pair, at least until the moment that it finds 
itself in state \(f>i,n— 2) and renewed relaxation may oc- 
cur. We see from this that spontaneous emission intro- 
duces a dephasing of the coherent oscillations of A/"(t) in 
different Landau-Zener pairs. This means that the neg- 
ativity computed for each pair is uncorrelated in phase 
from that of it's next neighbors. The total negativity is 
very sensitive to decoherence, the total negativity being 
smaller than the sum of the negativities over all Landau- 
Zener pairs. This example illustrates the subadditivity 
character of the negativity, in particular the requirement 
of phase correlation in superposition of entangled state 
parts. On this basis we may attempt to discuss in a 
broader sense the decay of entanglement seen in Fig. [2]d. 

As cooling ensues, a spread of the density matrix ele- 
ments occurs. This is evident from Fig. [8] which gives a 
matrix representation of the absolute values of elements 
of the density operator. Each square representing a vi- 
brational quantum number contains nine elements for 
electronic population and coherences. The result shown 
here corresponds to the conditions T = 6, gi =hgi = 1.34, 
ui = 0.03, A = 15 in units of 2tt MHz, and 77 = 0.1 
with n- m = 15. Note that here the negativity reaches a 
higher value, and does so faster than in Fig. [2b, where 



FIG. 8: The absolute density matrix elements are shown at 
a time t = 24 ps, when negativity maximizes for these pa- 
rameters, see text. The Landau-Zener model for negativity is 
confined to the elements at the corners of a single red square. 



.9i/52 = 10. Reasons are the two times higher two-photon 
Rabi- frequency and the increased coupling (21) which 
scales with 



Jn- m . Fig. ffl shows the density matrix at 
a time when the maximal negativity is reached (see in- 
set). In our Landau-Zener model the entanglement stems 
solely from the contributions of the elements of a sin- 
gle Landau-Zener pair, marked by the corners of the red 
square in Fig. [8| The negativity based on the trace norm 
takes into account all the density matrix elements from 
the master equation and we see from Fig. [8] that pop- 
ulation and coherences have already spread outside the 
initial Landau-Zener pair. 



C. Entanglement and spontaneous emission 

Returning to the full model simulation which gives us 
the results shown in Fig. [2] we evaluate the rate of emis- 
sion of spontaneous photons from the full density matrix 
using 



R(t) =r^p3n 3n (t) 



(33) 



Figure [9] shows the development of this rate during the 
cooling process for the three cases considered in Fig. [2] 
From Eq. ( |33[ ) we may estimate the total number of spon- 
taneously emitted photons (approximately 2 spontaneous 
photons per vibrational level cooled [12]) and also the 
time of emission. We may argue that at times r when 



Nphoton(T) = T / dt 

Jo 



P3n3n(t) 



(34) 
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climbs to the next integer value a spontaneous photon 
has left the trapped atom. The sequence of integer val- 

^pholon{ T ) 

1 2 345 10 15 20 



0.008 



To explore this situation at later times we show in 




0.000 



Fig. 11 the distribution of density matrix elements at a 
time when the negativity is near its lowest values. Again, 
red squares mark the elements of individual Landau- 
Zener pairs in the density matrix. It is evident that 
following the initial climb to a high negativity a drop 
into a region of low negativity ensues, which is caused by 
spontaneous emission and which reflects the incoherent 
distribution of density over many elements. Neverthe- 
less, substantial coherences appears within each square, 
as is evident by the strength of the elements at the off- 
diagonal corners. Incoherency is a consequence of the 
statistical character of spontaneous emission. Only it 
leads to a transfer of population from one square to it's 
lower neighbor and, hence, to cooling. 



FIG. 9: Rate of spontaneous emission during cooling for the 
three cases shown in Fig[2] Integral values of the accumulated 
rate (the number of spontaneous photons emitted) are given 
by the numbers on the top axis for r = 6 MHz and 200 kHz 
respectively. 

ues N photon (r) is marked along the top axis in Fig. [9] 
for the cases T = 6 and 0.2 MHz. We see that sponta- 
neous emission practically ceases when the cooling limit 
is reached. The connection between the oscillatory na- 
ture of the rate of emission of spontaneous photons and 
the rate of oscillations of the negativity in Fig. [2Jd can 
be understood by comparing these results directly, see 
Figure [10] The rate of spontaneous emission is perfectly 
out of phase with negativity and maximizes whenever the 
state \4>i,n—l) is most strongly populated, that is at ev- 
ery second minimum in negativity, in perfect agreement 
with the prediction of the Landau-Zener model. 



o - 

2 - 



I .0- 
12 - 
14 - 
16- 



0.1 






0.01 






r 


0.001 








■ I 4 








10 s 








18 - | Pinjm I 





0.5 




0.4 


o 


0.3 


II 






0.2 








0.1 




0.0 




10 r 




8 


| 






6 




4 






n = 


2 


u 






n 



r=200 kHz 
r=600 kHz 
1=6 MHz 




1000 





T=6MHz 


t 

7 ^\ 

! \ y 


T=600 kHz 




T=200 kHz 



2(H) 



400 600 
time (^s) 



800 



1000 



FIG. 10: The rate of spontaneous emission (lower figure) is 
out of phase with the oscillations in negativity (upper figure) 
and occurs preferentially at every second minimum of nega- 
tivity, when the state \<j)i,n — 1) is most strongly populated. 



FIG. 11: The absolute density matrix elements are shown 
at a time t = 400 fis, when negativity has dropped to near 
zero. Parameter as in Fig. [8] Populations and coherences 
have by now spread over a wide domain of the density matrix. 
The action of many neighboring Landau-Zener pairs is clearly 
evident as marked by the red squares. 



D. Rebirth of entanglement 

During the cooling process the atom is driven into the 
unique stationary state p s tat of the master equation pj). 
The final value Af(oc) of the negativity represents the 
entanglement of this state which has a dominant contri- 
bution from \(f>2, 0), the product of the dark state and the 
vibrational ground state of the trap with zero negativity. 
However, p s t a t contains a small amount of entanglement 
which is controlled by the off-diagonal transitions from 
the dark state. This can be seen from the dynamics of 
the fidelity [15] of the dark state, 



F(t) = (^2,01^)102,0), 



(35) 



9 



which is shown in Fig. 12 1. We find that the fidelity 
at the end of the cooling cycle (the stationary state) is 
approximately equal to 0.989 for Q\lg% = 10 and 0.964 
for 51/52 = 5 for the parameters used in Fig. [2] with 
r = 6 MHz. 




2 5 10 20 
t (ms) 



ftnlm P3n2m P3n3i 



FIG. 12: (a) The fidelity of the dark state 1 4>2, 0} as a function 
of time at F = 6 MHz for the parameters of Figs. [2] and 1 1 1 1 
respectively, (b) Matrix representation of the cooling limit 
(t = 6000 fj,s) for the conditions depicted in Figs. [8] and |11| 

The observed revival of the negativity at the end of the 
cooling process thus reflects the amount of entanglement 
in the final stationary state. To estimate the final nega- 
tivity for small values of T we determine the interacting 
ground state |xo) of the Hamiltonian H v to first order 
in the Lamb-Dicke parameter. Taking into account the 
coupling of the state \4>2,fy to the excited state |0i,l) 
(see Fig. [5]) through the matrix element 

^=<^,l|^Sl^,0)»-i^, (36) 

we find the ground state (normalized to first order in 77) 

A 



\Xo) 



102,0) 



2lj 



101,1) 



The negativity of this state is given by 
\rt\ w l-^l 9192 



(37) 



(38) 



In the limit T — !• this negativity should coincide with 
the value TV (00) for the stationary state negativity. In 
fact, this prediction is in quantitative agreement with the 
results for the full model in Fig. [2]d: For the parameters 
chosen in this figure Eq. (38) yields JV(oo) = 0.01 in 
very good agreement with the simulations results for T — 
200 and 600 kHz. For T=6 MHz the two-photon Rabi 
frequency is no longer much faster than the spontaneous 
rate and, hence, in this case the result of the numerical 
simulation in Fig. [2}d falls slightly below the prediction 
of Eq. p8| by about 20%. 



E. Entanglement control 

The discussion above shows that the entanglement in 
the trapped atom's degrees of freedom is imprinted by 



two-photon transitions between the two ground states. In 
the absence of an external radiation field the two ground 
states are however free from decay. We therefore expect 
that entanglement can permanently be imprinted in the 
trapped atom system. Here we explore this effect by 
switching off both lasers at a certain time during the 
cooling cycle. We do not consider the case that only 
one laser is turned off, as this will obviously lead to a 
complete repumping into one of the ground state levels 
by spontaneous emission with a steady-state negativity 
of zero. 

The turnoff of the laser fields starts at some time 
t ff and lasts for a time At, and is described by time- 
dependent Rabi frequencies gj(A(i)) = giX(t) with the 
following form for the function X(t), 



X(t) 



2 TT{t-t aii ) 

1 2At 



for t < t Q ff, 

for t ff <t <t ff- 
for t > t ff + At. 



-At, (39) 



To model a sudden switch off we choose At = 1/is. Two 
examples of simulation results for the full master equa- 
tion are shown in Fig. [l3j We notice that the negativity 
is practically frozen at its value at the turnoff time. A 
closer inspection reveals a minute decrease in negativity 
on the time scale of spontaneous emission, which is con- 
nected to the small contribution to the negativity from 
the excited electronic state. 




O.HI 
0.08 
' 0.06 
0.04 
0.02 
0.00 
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50 100 
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FIG. 13: A sudden turnoff of both lasers nearly preserves the 
negativity. The close-ups show that a minute fraction of neg- 
ativity suffers decay which stems from the tiny contribution 
of the excited electronic state. Parameters as in Fig. [2Jd. 



The observed freezing of negativity can be understood 
in general terms by means of the following argument. The 
Hamiltonian of the total system after switching off the 
lasers takes the form H cm + H c \, describing an uncoupled 
time-evolution of the internal and the external degrees of 
freedom. Since the contribution of the excited electronic 
state is very small under EIT conditions, the total time- 
evolution operator after switch-off is given by a tensor 
product of local unitary operators which leaves invariant 
any measure for entanglement. 

When turning off both lasers slowly a more adiabatic 
transition regime is observed, see Fig. |14| Now the neg- 
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FIG. 14: Adiabatic turnoff of both lasers over a time period 
of At = 50 /xs. The two examples shown refer to t s = 100 /is 
and t ff = 200 /us. 



ativity during turnoff mirrors th e de creasing Rabi fre- 
quencies in the coupling element (21), A 



This element controls the oscillations of negativity with 
a single Landau-Zener pair. The adiabatic decrease in 
the two-photon Rabi frequency diminishes this oscilla- 
tion, regardless of whether we are in the rising or de- 
creasing regime of negativity. Similar observations were 
made at T — 6 MHz, albeit at overall lower values of 
negativity. 



IV. DISCUSSION AND CONCLUSIONS 

We have discussed here in detail the entanglement dy- 
namics in EIT laser cooling of trapped atoms, employing 
the negativity as a measure for entanglement. We have 
demonstrated that essential features of the dynamical be- 
havior of the negativity can be explained in terms of a 
conceptually simple model which is based on Landau- 
Zener state pairs at avoided crossings in the spectrum of 
the full atom-field Hamiltonian. As we have seen by the 
comparison with numerical simulations of the full master 
equation, the model correctly describes the entanglement 
oscillations and their decay in the early stage of the cool- 
ing process, as well as the revival of entanglement towards 
the final stationary state of the atom. While these fea- 
tures can be explained qualitatively and quantitatively 
on the basis of the Landau-Zener model, a deeper under- 
standing of other features, like the observed power law 
decay of negativity and its precipitous drop in the tran- 
sient regime of cooling, is still lacking. This requires a 
more sophisticated treatment of the dynamical transfer 
of negativity between different Landau-Zener state pairs 
induced by spontaneous emission processes. In this con- 
text it may also be of interest to explore the situation 
when two or more initial vibrational states are driven 
into coherent oscillations of negativity. 

In Sec. MI El we have discussed a scheme for the freez- 
ing of the entanglement between internal and external 
atomic degrees of freedom by a rapid turnoff of the laser 



fields. This scheme may find further interesting applica- 
tions in the entanglement control of a much more general 
class of physical systems which are driven by external 
fields. Moreover, it seems worth studying in more detail 
other control methods involving an adiabatic turnon and 
turnoff, or suitably designed pulse shapes and sequences 
of pulses of the external driving fields. 

Finally, the theoretical picture developed here in order 
to understand the dynamical behavior of the negativity 
points towards fundamental distinctions between the en- 
tanglement which is generated by laser cooling of atoms 
in the EIT scheme on one hand, and by VSCPT on the 
other. We briefly discuss these distinctions in the follow- 
ing. 

According to the Landau-Zener model constructed in 
Sec. |IIB| the initial phase of the entanglement dynamics 
in EIT cooling can be understood as a periodic oscillation 
between the eigenstates 



X±) = 



1 

71 



(\<t>x,n-l)±\^ 2 ,n)) 



(40) 



of the trapped atom-laser field Hamiltonian H n at a 
Landau-Zener splitting. These states are maximally en- 
tangled and lead to a periodic oscillation of the negativity 
between zero and \ with a period T cnt — ir/AE which is 
determined by the energy gap AE at the avoided cross- 
ing. The periodic oscillation of the negativity is inter- 
rupted stochastically by spontaneous emissions. After a 
spontaneous emission event a build-up of entanglement 
in the next lower Landau-Zener pair ensues. As this con- 
tribution is however unrelated in phase to that of the par- 
ent pair, the subadditivity of the negativity leads to an 
overall decrease of the amount of entanglement as cooling 
proceeds. The dominant final state of the system which is 
reached in EIT cooling is the product state \<f>2, 0). How- 
ever, a small degree of entanglement appears in the final 
state of EIT cooling which is mainly due the coupling of 
the states \4>2, 0) and \<j>\, 1), as is indicated by the lowest 
diagonal blue arrow in Fig. [5] In the limit of T — >• a 
small negativity results from this coupling which is given 
by Eq. p|. 



For the VSCPT scheme of laser cooling the dynamical 
behavior of entanglement is markedly different. In this 
scheme a free atom is irradiated by counter-propagating 
orthogonally polarized laser fields. As long as real spon- 
taneous emission events are neglected, the motion of the 
atom is confined to the three-dimensional manifold F(p) 
which is spanned by the product states 



{|3,p), \l,p-k), \2,p + k)}. 



(41) 



where p denotes the atomic momentum, and k the pho- 
ton momentum. This fact is a consequence of the con- 
servation of the total momentum of the atom-laser field 
system, which does not hold for an atom trapped in an 
external potential. As a result of the finite widths of the 
dressed atomic states, the atomic wave function is rapidly 
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pumped into the nearly dark state 

|x(p)) = -^[|l,p-fc> + |2,p + fc>] (42) 

during the evolution within a particular manifold F(p)- 
This is a maximally entangled state of the internal and 
external degrees of freedom of the atom with a nega- 
tivity of h. The effect of spontaneous emissions is that 
the atomic state vector leaves the manifold F(p) to enter 
a new manifold F(p') in which it is then again driven 
rapidly to a state of the form ( 42 ) with the new mo- 
mentum p' . Within this scheme the cooling proceeds 



through the consecutive reduction of the atomic momen- 
tum to values smaller than that corresponding to the 
photon momentum, because the atoms are approaching 
a state which is completely decoupled in the limit p — > 0. 
In fact, in this limit the state 

|x(0)> = -^[|l,-fc> + |2,+fc>] (43) 



is reached which is a true dark state that cannot be re- 
excited by the laser fields. Thus we see that in the long 
time limit of VSCPT cooling the atom is nearly always 
in a state of maximal entanglement between the inter- 
nal and the external degrees of freedom. These states 
also exist after an instantaneous turnoff of the driving 
fields. The resulting atomic wavepacket spreads however 
in space at a rate proportional to the recoil momentum. 
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